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Abstract. A model for a solar flare, involving magnetic reconnection transferring 
flux and current between current-carrying magnetic loops connecting two pairs of 
footpoints, is generalized to include conservation of magnetic helicity during re- 
connection, as well as conservation of current at all four footpoints. For a set of 
force- free loops, with the ith loop having flux Fi and current h, the self and mutual 
helicities are proportional to the self and mutual inductances with the constant of 
proportionality determined by cu = Fi/fioli. In a constant-a model, the change 
in magnetic energy is proportional to the change in helicity, and conservation of 
helicity implies conservation of magnetic energy, so that a flare cannot occur. In 
a quadrupolar model, with ai > ai initially, ct\ increases and Q2 decreases when 
flux and current are transferred from loops 1 and 2 to loops 3 and 4. A model that 
conserves both current and helicity is constructed; it depends on the initial as, and 
otherwise is somewhat simpler than when helicity is neglected. 
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1. Introduction 

Solar flares are attributed to magnetic energy release in the solar corona. 
A generic flare model involves a pre-flare magnetic configuration relax- 
ing through magnetic reconnection to a post-flare configuration with a 
lower energy, with the magnetic energy difference identified as the en- 
ergy released in the flare. However, the data on pre- and post-flare mag- 
netic configurations are far from definitive in determining the details of 
the change in the magnetic configuration. As a consequence all detailed 
flare models rely on additional assumptions to constrain the pre- and 
post-flare configurations. One class of model involves reconnection be- 
tween two current-carrying magnetic flux loops (Machado et aL, 1988; 



Kundu et aL, 1995||Wang et aL, 1995l|Uchida et aL, 2003[ ). Such mod 



els are referred to as 'two-loop' models, although four loops connect- 
ing the two pairs of footpoints are actually involved; the alternative 
name 'quadrupolar' (Uc hida et aL, 20031 ) is used here. In such a model, 
magnetic reconnection involves transfer of magnetic flux and current 
between two pairs of loops. During a flare, the magnetic configuration in 
the corona changes, but there is insufficient time for the magnetic field 
below the photosphere to change. This is due to the change in magnetic 
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stress propagating at the Alfven speed, with the Alfven propagation 
time being shorter than the flare time in the corona but longer than 
the flare time for propagation to any significant depth below the photo- 
sphere. This leads to the important constraint that the flux and current 
passing through the photosphere at each footpoint is unchanged during 
a flare. This constraint rules out several otherwise possible models, 
including a global change from higher to lower a and annihilation of 
two oppositely directed currents. 

The non-potential component of the magnetic energy, due to the 
currents flowing in the corona, is described in terms of the currents 
and the self and mutual inductances. Current conservation implies 
that the changes that lead to magnetic energy release must be in the 
configuration of the current, rather than changes in the magnitude of 
the current. In a force- free corona, the current paths are along mag- 
netic field lines. A simplifying assumption that allows one to model the 
change in the current configuration is that there are fixed current paths 
(loops) in the corona, with a fixed set of self and mutual inductances 
describing these paths. The changes during a flare correspond to cur- 
rent being transferred between these paths ( |Melrose, 1997 ). The model 



involves a large number of parameters: ten (four self inductances, six 
mutual inductances) describing the geometry, four describing the initial 
currents in the four loops, and two describing the current and flux 
transferred. Simplifying assumptions are made to reduce the number of 
free parameters, allowing one to identify the most favorable configura- 
tions for a flare with the maximum magnetic energy difference between 
the initial and final configurations ( H ardy, Melrose and Hudson, 1 998 ) . 
Favorable configurations for magnetic energy release implied by this 
model are found to be correspond to favorable configurations found 
observationally for flares (Aschwanden at al, 1999). 



Magnetic helicity is conserved during reconnection (Taylor, 1986) 



and such conservation is believed to apply to solar flares 
Bcrgcr, 1999| |Priest, 199*9 ). Conservation of helicity implies that helic 
ity transported into the corona can only be rearranged in the corona, 
and must ultimately be transported out of the corona. Within the 
corona helicity can only be transferred from one set of field lines to an- 



evstov, Canfield and Zirin, 1 996; 



other, from one scale to another or between twist and writhe ( Yousel and Brand enburg, 2 003 ) 
The transport of helicity into and out of the corona has received con- 
siderable attention recently, both for injection through the photosphere 



(Bergcr and Ruzmikin, 2000; 


Kusano et al, 2002||Demoulin and Berger, 2003| 


Magara and Longcope, 2003 ; 


Welsch and Longcope, 2003||Zhang and Low, 2003 ] 



through the photosphere ( |van Ballegooijen and Ma rtens, 19891 ). The 
time for transport of helicity out of the corona is long compared with 
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the reconnection timescale, and hence helicity should be conserved 
on the flare timescale. Consequently, release of energy during a flare 
requires that the magnetic configuration change from one of higher 
magnetic energy to one of lower magnetic energy, with the net helicity 
of both configurations being the same. 

In this paper conservation of magnetic helicity is included in the 
reconnecting loop model. This provides an additional constraint on the 
flare model: magnetic flux and current at the footpoints, and helicity 
in the corona are all assumed to be conserved between the pre- and 
post-flare states. Note that conservation of helicity in a flare does not 
necessarily imply that the helicity is conserved between the observed 
pre- and post-flare configurations in the corona. For example, in erup- 
tive flares the helicity carried off by CMEs ( Rust, 1999 ) needs to be 
included in the helicity budget. 

The requirement of helicity conservation is included in a force-free 
multi-circuit model in section 2. The quadrupolar flare model is summa- 
rized in section 3. In section 4, implications of conservation of helicity 
are considered in a model for differential flux and current transfer are 
considered. The model is generalized to an arbitrary transfer of flux 
and energy in section 5. The results are summarized and discussed in 
section 6. 



2. Multiple circuit approximation 

The magnetic energy may be written in terms of the current. One starts 
from the general form 

£mag = ^|d 3 xJ(x)-A(x), (1) 

where J(x) is the current density, and writes the vector potential (in the 
Coulomb gauge) in terms of the current by solving Poisson's equation: 

A(x) = MQ /dV^-. (2) 

The resulting double integral, 

^ = a/^JWiM (3) 

2 J |x — x'l 

is the basis for a multiple circuit model (e.g., Jackson 1975). The current 
configuration is approximated by a set of discrete currents, I{ with 
i = 1, 2, . . ., and the integral © is approximated by 

£mag = £ £ M^ld, = lE L ^ + E M iohIj, (4) 

ij i i<j 
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where Li = Ma is the self-inductance for the ith current and M« = M„ 
for i ^ j is the mutual inductance between the ith and jth currents. 
By choosing a sufficiently large number of discrete currents, the mul- 
tiple circuit model may be used to obtain approximations of arbitrary 
accuracy. However, the model is useful in the present context only if 
the configuration of interest can be approximated by a small number 
of discrete currents. 

The magnetic helicity, H, is defined by 

H = J d 3 xB(x) • A(x), (5) 

and it may be expressed in terms of the current using © twice, with 
B(x) = curl A (x). A simplifying assumption made here is that the only 
contribution of interest is for force-free configurations of the magnetic 
field. This corresponds to 

J(x) = ^Mb(x), (6) 
Mo 

with a(x) constant along each current line, which coincides with the 
magnetic field line. Then the integral © with © is of the same form 
as the energy integral (^Q) . In the multiple circuit model the relation @ 
implies a relation 

U = -Fi (7) 
Mo 

between the current and the magnetic flux associated with the ith loop, 
and where on is an appropriate average of a(x) over the cross-section 
of the loop. Then the analogy with the magnetic energy implies that 
the helicity may be written in the form 

# = E m ^ J i = E + E 2m « 7 i 7 J . ( 8 ) 

ij i i<j 

with the term involving li = mu describing the self helicity and the 
terms involving rriij = rriji describing the mutual helicity. These terms 
are related to the self inductances Li and the mutual inductances M^j 
by 

H — j >'Hj — > "ti — , > \&1 

oti aij (\, + n ; 

where the symmetry property rriij = m ji is imposed without loss of 
generality. 

Alternatively, the expression (JSJ) for H may be expressed in terms 
of the fluxes using (jJJ). This gives 

H = J2 C i3 F i F J = E T i F i + E 2 ^3 F i F 3 . ( 10 ) 

ij i i<j 
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with the term Tj = Ca describing the self helicity and the terms in- 
volving dj = Cji describing the mutual helicity. The form ()10j) for H 
was written down by Berger (1998; 1999). The coefficients in (|1U|) are 
related to the self inductances Lj and the mutual inductances My by 

, Oiij _ Gin ~T~ Gi i* _ _ , . 

Ti = — L h dj = 3 - Mij. 11 

A*0 IM) 

The quantity Tj may be interpreted as the number of twists (angle of 
the twist divided by 2tt) along the ith loop. 

Although explicit expressions for the coefficients Tj and dj have 
been presented in the literature, these are not satisfactory for present 
purposes. There are several properties that an acceptable form of Cij 
should satisfy: (i) dj = £ji should be symmetric in interchange of the 
loops (any asymmetric part has no physical significance), (ii) dj should 
have its maximum positive value when the loops are identical, parallel 
and coincident, when dj should be equal to Tj = Tj (the currents and 
fluxes of the two loops then add, (iii) dj should have its maximum 
negative value when the loops are identical, antiparallel and coincident, 
when dj should be equal to — Tj = —Tj (the combination of the two 
loops is null) , and (iv) dj = when the two loops are orthogonal. The 
available calculations of the mutual helicity involve first showing that 
the rate of change of the mutual helicity between loops (flux tubes) i 
and j is 

dH V _ 2 r> FF ( U ) 

where Ojj is 'the angular velocity of the jth flux tube about the ith' 
(Welsch and Longcopc, 2003), and then integrating from some specific 



motion starting from an initial condition where Hij is zero (Berger, 1998 ) 



The resulting expression (Berger, 1998 ) appears to satisfy none of the 



requirements (i)— (iv). This difficulty is implicit in the lack of symmetry 
between i, j in the interpretation of Qij in Q12|). and in the lack of 
symmetry in the manner in which the integral over time is performed. In 
a detailed theory, in which an explicit expression for the mutual helicity 
is required, this difficulty must be addressed explicitly. However, the 
arguments in the present paper do not rely on any explicit expression 
for Cij, and this difficulty is not discussed further here as it does not 
affect any of the conclusions. 

Nevertheless it is appropriate to identify an explicit expression for 
the mutual helicity that satisfies the requirements (i)-(iv). A semi- 
empirical result follows by modifying an approximate expression for 



the mutual inductance (Melrose, 1997), which depends only on the 
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Figure 1. The quadmpolar model: (a) the two initial loops are shaded lightly and 
the two final loops are dark; (b) in projection on the photosphere, the initial loops 
are the solid lines, the final loops are the dotted lines, and the faint curves denote 
transferring flux and current. [After Melrose (1997)] 



geometry of the current paths. This semi-empirical result is 

/ 4 \ 3/2 

^ = (W /2 ™» tl [ (ai + ^ +(tij ) ■ (13) 

where 8ij is the angle between loops i and j, such that 6ij = when 
the loops are parallel, 0y = tt when they are antiparallel, and 0« = tt/2 
when they are orthogonal, and where dj , a,- are the lengths of the loops 
(the major radii in a toroidal model) and dy is the distance between 
their centers. It is desirable that an expression for dj that includes the 
important features incorporated into be derived rigorously. 



3. A quadrupolar flare model without helicity conservation 



The flare model adopted here (Melrose, 1997), cf. figure 1, is based 



on the change in the magnetic energy in the form (py) as the result of 
a change in the magnetic configuration during a flare. An important 
assumption in the model is that only the initial (pre-flare) and final 
(post-flare) states are of direct interest, and in both cases the energy 
may be written in terms of the currents and inductances. (The induc- 
tances are changing during the reconnection process, and an energy 
integral then does not exist in general.) 
The energy released is identified as 

Eq = -Emag — -Erriagi -^mag = 2 ^ij^i^ji (14) 

ij 

where the unprimed and primed quantities denote the pre- and post- 
flare configurations, respectively. The assumption that the geometry 
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does not change corresponds to = Mij. A flare is possible only if 
there is free energy to drive it, and this requires Er > 0. In principle, 
(|14|) could be used to calculate the energy release as a result of any 
change in the magnetic structure and current distribution in the corona, 
simply by modeling the pre- and post-flare fields with as many loops 
as desired. 

A generic model involves transfer of flux and current between pairs 
of footpoints. Let the ith pair of footpoints be denoted i±. A partic- 
ular loop connects two footpoints, and let the loops be denoted by 
the footpoints that they connect. Thus loop connects footpoints 
i + and j_, and has a magnetic flux Fi + j_ and a current Both 
flux and current must be conserved at each footpoint in such a model. 
Conservation of current at the footpoint i+, say, requires 

3- 3- 
and conservation of magnetic flux requires 

Fi + =F' i+ , F i+ =Y,F i+j _, ( 16 ) 

3- 3- 

where the sums are over all the footpoints j_ connected to £+ by a 
loop. There is one such constraint for each footpoint. 

In practice, we have no direct information on the reconfiguration of 
the magnetic field, and we need to make major simplifying assumptions 
to formulate a flare model. The assumptions made here are: 

(a) The current paths (loops) are predetermined, and the changes 
involve redistributing the current between these paths. In other 
words, it is assumed that the transfer of flux and current is between 
loops with a fixed geometry and this geometry (which determines 
the inductances) does not change as a consequence of this transfer, 
so that one has M'- = My. 

(b) A flare is modeled in terms of current and flux transfer between 
two pairs of footpoints, so that only four loops are involved. 

Note that in the quadrupolar model there are four loops, thought of 
as two initial loops and two additional final loops formed by making 
new magnetic connections between the four footpoints. Let loops 1 + 1_. 
2 + 2_, l+2_, 2 + l_ be labeled more simply as loops 1, 2, 3, 4, and let the 
initial and final currents flowing in these loops be denoted by unprimed 
and primed quantities, respectively. It is assumed that a current AI 
and a flux AF are transferred from the initial loops (1 and 2) to the 
final loops (3 and 4). Conservation of current implies I[ = I\ — AI, 
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I' 2 = I 2 — AI, 13 = 73 + AI, 1^ = Z4 + AI, and conservation of flux 
implies an analogous set of relations. It is convenient to write the energy 
released in the flare as (Melrose, 1997) 

E R = RAI + M IR (A/) 2 , (17) 

with the coefficients given by 

R = M\ CS I[ + M 2 LCS ^ - M 3 LCS 7 3 - M4 cs / 4 , 

Mf cs = M jx + M j2 - M j3 - M jA , (18) 
M IR = UL X +L 2 -L 3 - L a ) + M 12 - M 34 , (19) 



where LCS labels 'like current separation' terms and IR labels 'irre- 
ducible reconnection' terms. 



4. Infinitesimal change conserving current and helicity 

It is helpful to consider the case of an inifinitesimal reconnection (AI —* 
dl, AF —* dF) that conserves both current and helicity. This implies 
a change in the helicity due to a change in the currents, ij, rather than 
to a change in the geometry of the current (or magnetic field) pattern, 
as assumed in (|12|). 

4.1. Change in the magnetic energy 

An infinitesimal change Ii — ► 1$ + dli in the currents implies an in- 
finitesimal change in the magnetic energy: 

dE m ^ = Y,M i] I l dI ] . (20) 

ij 

Conservation of current in the quadrupolar model requires dl\ = dl 2 = 
—dl 3 = — dli = —dF The flux also changes, with dF\ = dF 2 = —dF 3 = 
-cLF 4 = -dF. The becomes 

d^mag = - ^(Mali + M i2 I 2 - M i3 I 3 - Mih)dF (21) 

i 

One requires dE ma _ g < for magnetic energy to be released to drive a 
flare. 
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4.2. Change in the as 
Suppose one writes 



Equal amounts of current and flux are transferred from the two initial 
loops to the final loops, and if the initial loops have different as, then 
one expects the value of ao to be a weighted average of ai and ot2- 
Moreover, if the two final loops are created by the reconnection then 
one expects them to have 03 = = ao; this is because the same ratio 
(ao/Vo) of current to flux is transferred from both initial loops to both 
final loops. 

Assuming a% > ao > 0,2, a notable consequence is that a±, ct2 do 
not remain constant. Writing /^(I, + dZj) = (a^ + dcti){Fi + dFi), one 
finds 



In particular, one has dot\ = — ai(l — a\/ oto)dl / 1\ > 0, and da2 = 
—02(1 — ctz/aQjdl /I2 < 0. The loop with the larger (smaller) a has its 
a increased (decreased). Superficially it might appear that this result 
violates Taylor's theorem, which implies that the system should relax 
to one of constant a. In fact, the overall mean value of a does decrease 
as a result of the current and flux transfer. 

4.3. Infinitesimal change in the helicity 

The change in the as implies that the coefficients my in the expression 
Q for the helicity change. Using (|9"|). with dMij = by hypothesis, the 
change in the helicity becomes 



where ()23j) is used. 

The change in the energy Q2U|) is subject to the condition dH = 
with dH given by (|24|). Consider the identity dE milg = dE majg + 
CdH, which is satisfied for all constants C. Choosing C such that the 
condition dE mSug = is manifestly satisfied for constant-a, that is for 
aj = ao for all i, corresponds to C = —ao/2fiQ. Then in place of (|2U|) 
one has 




(23) 




(24) 




(25) 
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Assuming 03 = = ao, in place of (|2*Tj) one has 

dEu* = ~ £ ( ^f^M ah + ^p°M l2 I 2 ) dl. (26) 

Note that one can rewrite (J2H1) m term of the fluxes: I\/ct\ = F\/hq, 
I 2 /a 2 = Fz/fio, dl/ato = dF/n . 



5. Flare model that conserves helicity 

The foregoing discussion of the differential case suggests how to proceed 
in the more general case where there is an arbitrary transfer of current, 
AI, and flux, AF = fi AI/a . 

5.1. Generalization of the quadrupolar model 

Provided thet helicity is conserved one may introduce an equivalent 
energy by adding an arbitrary constant, C, times the helicity to the 
energy: then the change in the equivalent energy is equal to the change 
in the energy. As in the differential case it is appropriate to choose the 
constant C = — ao/2^o m defining the equivalent energy: 

Fraag = F ma g — H, Ef[ = E ma g E m , (27) 

where conservation of helicity implies that the additional term oc H 
does not affect the value of the energy released. The choice (|27|) is 
made to ensure that the requirement E$ = in a constant-a model be 
manifestly satisfied. 

The change in the energy at constant helicity and constant current 
then reduces to 

^fl = -r{ ( l ~ — ) t 2L i + M i2 - M i3 ~ M u]h 



4 I V ai 

«2 



00 

Q l 



[2L 2 


+ Mia 


-M 23 


- M 24 ]I 2 


[2L 3 


+ M 34 


-M13 


- M 23 ]I 3 


[2L 4 


+ M34 




- M 24 ]h 



(28) 



One may compare (|28[). which incorporates helicity conservation, 
and ()17|1 with ()19() . which does not include helicity conservation. For 
the purpose of the comparison one needs to make the replacements 
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I[ -> h - AI, I' 2 -> I 2 - AI in (dZJ) with (HU). Inspection then shows 
that ()28j) is somewhat the simpler of the two expressions, in particular 
because there in no term oc (AI) 2 . The way the helicity constraint 
is used, to require that the energy change be zero in the constant-a 
case, has the associated consequence of removing terms oc (AI) 2 . In 
this sense, inclusion of the helicity constraint simplifies the model by 
relating the final configuration to the initial configuration, allowing one 
to eliminate explicit dependence on the final configuration. 

5.2. Model involving creation of final loops 

The simplification introduced by imposing conservation off helicity be- 
comes more apparent when one considers the simpler case in which 
the two final loops are created by the reconnection. This model implies 
13 = 14 = 0, I' 3 = I4 = AI, and it also implies 03 = «4 = «o- Hence, 
the terms oc 13,14 drop out for two separate reasons (13,1a, — > and 
03 = <xi — > cto). This suggests that the dependence of the result on 
the final configuration is now almost completely implicit in the helicity 
constraint. 

It is further plausible to assume that the characteristic Qo of the 
flux transferred is the mean of the initial as: 



The most favorable configuration for a flare is to be identified with the 
configurations that maximize Eq. 

The sign of Er is determined by the sign of a\ — a<i and the sign 
of the quantity in curly brackets in ()3U|) , For the terms involving the 
self inductances, which are the largest terms in many cases, to give 
£"fi > requires L\I\ja\ > L2/2/«2- Hence, for a\ > 02 the favorable 
requirement for a large energy release is -L1/1 3> L2I2, that is, the 
dominant initial loop (L\I\ 2> L2I2) needs to have the larger twist per 
unit length (ct\ > a^)- 

A simple model corresponds to all the current and flux in loop 2 
being transferred to the final loops. This corresponds to I2 = AI and 
a?2 = cto. Neglecting the term oc I2I0L2 in (|3U|) . and maximizing the 
resulting expression at fixed F2 = AF implies «2 = g^i- Then (|3UI) 



ao = \(a\ + 02). 
With these assumptions (f2"5|) simplifies to 



(29) 




(30) 
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simplifies to 

E & = i- [2Li + Mia - M 13 - M 14 ] h AI. (31) 
lb 

Ignoring the mutual helicities, the change in the magnetic energy in 
loop 1 is L±IiAI, and 1)3 lj) implies that one eighth of this energy 
becomes available as free energy to drive the flare. Comparison of this 
result with the corresponding result for this simple model for the case 
without the helicity constraint, one finds that (|17|) with 1)19)1 allows a 
greater fraction of the energy to be released as free energy. Hence, the 
helicity constraint precludes certain final magnetic configurations that 
might otherwise appear favorable for maximizing the energy release. 

6. Discussion and conclusions 

The primary objective of this paper is to include conservation of he- 
licity in a quadrupolar flare model that is based on conservation of 
current and flux at two pairs of footpoints. To achieve this objective, 
the first step is to write the helicity in a form analogous to the energy, 
which is described by a collection of loop currents and their self and 
mutual inductances. The corresponding self and mutual helicities are 
proportional to the inductances, with the constants of proportionality 
involving the as associated with the currents. 

Several implications of conservation of helicity are identified. First, 
in a constant a model, conservation of helicity implies conservation of 
magnetic energy, and hence a flare is impossible, subject to the proviso 
that the fluxes and currents at the four footpoints do not change. 
Second, conservation of helicity precludes a model in which transfer of 
current and flux from the two initial loops through reconnection creates 
two new final loops. Effectively, the final loops must already be present 
before current and flux can be transferred to them. How this difficulty 
is to be overcome in the model is not entirely clear. One suggestion is 
that the two new loops form by splitting the initial loops with the split 
parts moving apart into their final positions. Third, the expression of 
conservation of helicity may be used to rewrite the change in magnetic 
energy in a manner that explicitly satisfies the requirement that the 
change in energy be zero when all the as are the same. This results in an 
expression of the change in magnetic energy which is proportional to the 
difference in the initial as, but otherwise has a relatively simple form. 
In effect the constraint on the allowed final configurations imposed 
by helicity conservation allows one to write the change in a way that 
superficially does not depend on the final configuration. Fourth, a flare 
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is allowed only if the stronger initial loop has the larger a, and the a of 
the post-flare (reduced flux and current) loop is higher than the initial 
value. 

The investigation in this paper is far from a complete discussion of 
the implications of conservation of helicity on the quadrupolar model. 
Three specific points clearly need further investigation. First, the appli- 
cation of the quadrupolar model to the interpretation of data on specific 
flares ( Aschwanden at ai, 1999 1 needs to be reconsidered, particularly 
in view of the importance of the ratio of the as in the two initial loops 
in the helicity-conserving form of the model. Second, a problem iden- 
tified in the course of this investigation concerns explicit approximate 
expressions for the mutual helicity. Available expressions (Berger 1998, 
1999; Priest 1999) do not satisfy four plausible requirements, (i)-(iv), 
discussed in section 2. A semi-empirical formula that satisfies these 
requirements is written down in (|13j) . Direct calculation of the mutual 
helicity needs to be reconsidered. Third, the model does not include 
magnetic structures that are unconnected to the photosphere, and such 
structures may form as a result of reconnection, and be identified with 
CMEs. 
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